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Abstract. While the need for accurate calculation of nucleosynthesis and 
the resulting rate of thermonuclear energy release within hydrodynamic mod- 
els of stars and supernovae is clear, the computational expense of these nu- 
cleosynthesis calculations often force a compromise in accuracy to reduce the 
computational cost. To redress this trade-off of accuracy for speed, we present 
an improved nuclear network which takes advantage of quasi-equilibrium in 
order to reduce the number of independent nuclei, and hence the computa- 
tional cost of nucleosynthesis, without significant reduction in accuracy. In 
this paper we will discuss the first application of this method, the further 
reduction in size of the minimal a network. The resultant QSE-reduced a 
network is twice as fast as the conventional a network it replaces and requires 
the tracking of half as many abundance variables, while accurately estimating 
the rate of energy generation. Such reduction in cost is particularly necessary 
for future generation of multi-dimensional models for supernovae. 


1. Introduction 

Examination of the process of silicon burning reveals that the nuclear evolution is dom- 
inated by large groups of nuclei in mutual equilibrium. During silicon burning, these 
quasi-equilibrium (QSE) groups, first discussed by Bodansky et al (1968), form well 


in advance of the global Nuclear Statistical Equilibrium (NSE) The QSE groups also 
persist after hydrodynamic changes have fragmented NSE (see Hix & Thielemann 1998, 
Meyer et al 1998). The existence of these groups greatly reduces the number of indepen- 
dent abundances, suggesting a way to reduce the computational cost of silicon burning 
without significant loss of accuracy 

Tracking the nuclear evolution from helium burning through to NSE requires a 
network that includes nuclei from a-particles to Zn. Silicon burning presents a particular 
problem as material proceeds from silicon to the iron peak nuclei not via heavy ion 
captures but via a chain of photodisintegrations and light particle captures. The minimal 
nuclear set which can follow this evolution is the set of a-particle nuclei; a, 12 C, 16 0, 
20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 40 Ca, 44 Ti, 52 Fe, 56 Ni, 60 Zn. For convenience we will name this 
set T and refer to its abundance as Y* . Silicon burning in fact presents a larger problem 
as the nuclear flow from silicon to the iron peak nuclei does not proceed through nuclei 
with N=Z, especially when significant neutronization has occurred (Hix & Thielemann 
1996). However, for nucleosynthesis calculations within hydrodynamic models, such 
compromise is often made necessary by the computational limitations. Furthermore, the 
small size of the a network (14 nuclei and 17 reactions) makes application of QSE to 
a-chain nucleosynthesis a pedagogically useful example 


2. Network Basics 


From a set of nuclear abundances, Y (= n/piV^) and the rates for the reactions which link 
them, one can calculate the time derivatives of the abundances, Y. With these deriva- 
tives, the abundances of the included nuclei are evolved, at time t, and over timestep 
At, according to the implicit prescription 


Y(t + At)-Y(t) 
At 


= Y{t + At) 


Solving Eq. 1 is equivalent to finding the zeros of the set of equations 
Z(t + At) = ^ + Z1M _ Y{t + At) = 0 . 


(1) 

( 2 ) 


This is done using the Newton-Raphson method (see, e.g Press, et al 1992), which is 
based on the Taylor series expansion of Z(t + At), with the trial change in abundances 
given by 

dZ{t + At) 


-l 


AY = 


3Y{t + At) t 

where dZ/dY is the Jacobian of Z. Iteration continues until Y(t + At) converges. 


( 3 ) 


3. The QSE-reduced a Network 

The objective of the QSE-reduced a network is to evolve Y 7 (and calculate the resulting 
energy generation) in a more efficient way. Under conditions where QSE applies, the 



existence of the silicon and iron peak QSE groups allows calculation of these 14 abun- 
dances from 7. For the members of the silicon group ( 28 Si, 32 S, 36 Ar, 40 Ca, 44 Ti) and the 
iron peak group ( 48 Cr, 52 Fe, 56 Ni, 60 Zn) the individual abundances can be calculated by 


Yqse,Si( A Z ) 


C{ A Z ) 
C( 28 Si) 


A — 28 

y( 28 Si)y~ 


and Yqse, Ni( a Z) 


C( A Z) 
C( 56 Ni) 


Y( 56 Ni)Ya 4 , (4) 


where we have defined 


C( A Z) 


G( A Z) / pN A \ 

2 A { 9 ) 


A-l 


A 5 exp 



and 


0 = 


( m u k B T \ V 2 

V 2nh 2 ) 


( 5 ) 


G( A Z) and B( A Z) are the partition function and binding energy of the nucleus A Z, Na 
is Avagadro’s number, k B is Boltzmann’s constant, and p and T are the density and 
temperature of the plasma. (A — 28)/4 and ( A — 56)/4 are the numbers of a-particles 
needed to construct A Z from 28 Si and 56 Ni, respectively. 

Thus, where QSE applies, Y 7 is a function of Y 71 , where the reduced nuclear set 
7 Z is defined as a, 12 C, l6 0, 20 Ne, 24 Mg, 28 Si, 56 Ni, and we need only evolve Y n . It 
should be noted that Woosley et al (1973) and Hix & Thielemann (1996) have shown 
that 24 Mg is ordinarily a member of the silicon QSE group, but for easier integration 
with a conventional nuclear network, we evolve 24 Mg independently. 

While Y n is a convenient set of abundances for calculating Y r , it is not the most 
efficient set to evolve, primarily because of the non-linear dependence on Y a . Instead we 
define Y G = [Y aG , Y( 12 C), Y( 16 0), Y( 20 Ne), Y( 24 Mg), Y StG , Y FeG ] where 


YaC 


YsiG 

Yp eG 


Y a + E E 

:€Si group 16 Fe group 

i€Si group 

E r.. 

tgFe group 


( 6 ) 


Physically, Y qG represents the sum of the abundances of free o-particles and those a- 

particles required to build the members of the QSE groups from 28 Si or 56 Ni, while Ys lG 

and YpeG represent the total abundances of the silicon and iron peak QSE groups. 

Corresponding to this reduced set of abundances Q is a reduced set of reactions, 

with quasi-equilibrium allowing one to ignore the reactions among the members of the 

QSE groups. Unfortunately, the rates of these reactions are functions of Y^, and are not 

easily expressed in terms of Y g . Thus, for each Y G , one must solve for Y n and, by Eq. 4, 
—* _ —* —* ~ 

Y , in order to calculate Y G which is needed to evolve Y G via Eq. 1. Furthermore, Eq. 3 
requires the calculation of the Jacobian of Z, which can not be calculated directly since 
Y s can not be expressed in terms of Y g . Instead we find it sufficient to use the chain 
rule, 

d fg _ qy g dY n 
dY G ~ dY *■ d? G 

to calculate the Jacobian. Analytically, the first term of this chain rule product is easily 
calculated from the sums of reaction terms, while the second term requires implicit 
differentiation using Eq. 6. 
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Figure 1. Temporal evolution of ihe group mass fractions for T» = 5 and 
p, = 10 9 gcm~ 3 . 

4. Explosive Silicon burning with the QSE-reduced a network 

Because the products of hydrostatic silicon burning are trapped deep in the potential 

well of their parent star, it is only by explosion that the interstellar medium is enriched 

in intermediate mass and iron peak elements. Thus the ability to model explosive silicon 
burning is of particular importance, as even material which reaches NSE undergoes 
significant changes as it cools. For brevity sake, we will here concentrate on explosive 
burning, and refer the reader to Hix, Khokhlov, Wheeler & Thielemann (1998) for a 
more complete discussion of the application of this method. 

To model silicon burning occurring as a result of shock heating, we will follow the 
approximation introduced by Fowler Sc Hoyle (1964). Therein a mass zone is instanta- 
neously heated by a passing shock to some peak initial temperature, Tg,, and density, p,, 
and then expands and cools adiabatically, with the timescale for the expansion given by 
the free fall timescale, thd = (24nGp)~^ 2 = 446pg '^ras. 

Figure 1 follows the evolution of the group mass fractions for an example of this 
explosive burning model starting from Tg, = 5 and p, = 10 4 * * * * 9 gcm -3 . Adiabatic cooling, 
which drops Tg below 4 after 9.4 ms and below 3 after 21.6 ms, freezes out the nuclear 
reactions before NSE is reached, resulting in incomplete silicon burning , as discussed by 
Woosley et al (1973). The ev'olution of the mass fractions by the QSE-reduced network 
matches well with those evolved by its conventional counterpart, until Tg approaches 3. 
Comparison of columns 2 Sc 3 of Table 1 reveals that the individual abundances also 
agree quite well With Tg— 4, the individual abundances, even those as small as 10 — 9 , 
agree to within 5%. 



Table 1. Comparison of network abundances for Tg t =5 0 and p,= 10 9 gcm 3 . 


Time (ms) 

t 9 

8.77 

4.07 

17 7 
3.29 

255 

0.01 

Nucleus 

Y net 

Yqst 


Y 

1 qse 

Ynet 

“He 

7.90xl0“ 7 

7.82xl0~ 7 

1.04x10-® 

1.01 xlO" 8 

1.94xl0 -14 

12 C 

1.96xl0 -7 

1.96xl0- 7 

3.99x10-® 

3.23x10-® 

3.90x10-® 

16 0 

7.34 xlO" 7 

7.39 xlO" 7 

5 26xl0 -7 

5.07xl0" 7 

5 27xl0~ 7 

20 Ne 

2.63 xlO" 9 

2.63xl0" 9 

1.69x 10 -10 

1.48xlO- 10 

9 88X10 -11 

24 Mg 

2 24x10“® 

2.26x10-® 

2.88xl0“ 7 

2.98xl0- 7 

2.80xl0" 7 

28 Si 

7 65 xlO" 3 

7.76 xlO" 3 

7.58xl0 -3 

7.86xl0 -3 

7.58xl0~ 3 

32 s 

4.93xl0“ 3 

4.96X10- 3 

5.16xl0 -3 

5.15xl0- 3 

5 16xl0 -3 

36 Ar 

1.43xl0- 3 

1 42 xlO -3 

1.27xl0 -3 

1.21xl0- 3 

1.27X10" 3 

40 Ca 

132 xlO" 3 

130 xlO" 3 

1.32xl0 -3 

1.22xl0" 3 

1.32xl0 -3 

44 Ti 

7.07x10-® 

6 90x10“® 

1.96x10-® 

1.72x10-® 

1.69x10-® 

48 Cr 

5 89x10-® 

5.58x10-® 

4.40x10-® 

1.19x10-® 

4.40x10"® 

52 Fe 

7.17xl0 -4 

6 93 xl0“ 4 

6 33xl0 -4 

3.05xl0~ 4 

6 33 xlO- 4 

56 Ni 

8.63 xlO" 3 

8.60 xl0~ 3 

8.73 xl0“ 3 

9.04xl0- 3 

8.73x10-® 

60 Zn 

6.18x10-® 

6 06 xlO" 8 

3.26x10-® 

3.23x10-® 

4.38xlO- 10 


As the matter continues to cool, many of the photodisintegrations responsible for 
QSE freezeout between Tg = 3.5 and 3, fragmenting the large QSE groups into smaller 
grouplets and individual nuclei as the remaining free a-particles are captured. This 
reduces the viability of a QSE based approach. While at T 9 =3.3, the group abundances 
of the silicon and iron peak groups (which account for 99.9% of the mass) calculated by 
the QSE-reduced a network differ by less than 1% from the freezeout abundances of the 
conventional network, individual abundances are beginning to show larger variations. 
Comparison of columns 4 & 5 of Table 1 reveals that these variations are largest near 
the group boundaries. However, among the dominant nuclei the abundances are good 
to ~ 10%. Comparison of the full network abundances at T 9 =3.3 (column 4) with 
those at freezeout (column 6) reveals that only the free a-particle abundance, and those 
abundances comparable in size, are strongly affected by these final a-particle captures. 
Since the post-freezeout captures do not strongly affect the abundances, the abundances 
of the QSE-reduced network, frozen when Tg ~ 3.5, provide reasonable estimates of the 
final abundances. 

Figure 2 displays the evolution of another example of this explosive silicon burning 
model with T 9l = 6 and p = 10 7 gem -3 . The lower density in this case results in a slower 
expansion, with Tg reaching 5, 4, and 3 after 77, 172, and 293 milliseconds, respectively. 
Though the rate of cooling is relatively slow, the temperature drops too quickly for the 
large a-particle abundance to be completely incorporated into heavier nuclei, resulting 
in the a-nch freezeout of Woosley et al (1973). Even with this large over abundance 
of a-particles, the QSE-reduced network reliably tracks the group abundances until Tg 
approaches 3. 

As the material cools, the large a-particle abundance and resulting flow upward 
into the silicon group prevents the abundances within the silicon group from declining 
as rapidly as QSE requires, disrupting the QSE groups. Comparison of abundances cal- 
culated by the conventional a-particle network, at T 9 =4, with those of its QSE-reduced 
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Figure 2. Temporal evolution of the group mass fractions for r 9 , = 6 and 
p, = 10 7 gem -3 . 


counterpart (columns 2 & 3 of Table 2) shows the beginning of this process, with the 
abundances of 28 Si and 32 S much larger than QSE would predict. In spite of this break- 
down in QSE, the abundances predicted for the more abundant members of the group 
agree well. As temperature continues to drop, the disruption of the QSE groups by the 
large a abundance continues. By the time Tg drops to 3 (columns 4 & 5 of Table 2), 
the under prediction by QSE of the smaller group abundances also affects the iron peak 
group. However the abundances of the dominant nuclei, and hence the energy produc- 
tion, are still in good agreement. As the remaining photodisintegrations freezeout, the 
continued capture of the large abundance of o-particles results in significant abundances 
changes (column 6 of Table 2). However, even the conversion of 1.5% of the 56 Ni abun- 
dance to 60 Zn over several hundred milliseconds does not significantly affect the energy 
release of silicon burning. 


5. Conclusion 

We have shown that a QSE-reduced a network can be used as a replacement for full 14 
element a network when modeling silicon burning without significant errors in energy 
generation or nucleosynthesis. Such an approach offers a significant reduction in the 
computational time spent on nucleosynthesis, a factor of 2 for this small system of 
equations. It also offers a reduction in the number of nuclear variables which must 
be included within a hydrodynamic model, a matter of concern for models with large 
numbers of zones or grid points. We plan to apply this method of QSE-reduction to 





Table 2. Comparison of network abundances for Tg,= 6 0 and p,=10 7 gcm 3 . 


Time (s) 

t 9 

.174 
3 98 

.290 

3.02 

2.67 

0.01 

Nucleus 

Ynet 

Y 

1 qse 

Ynet 

Yqse 

Y ne t 

4 He 

1 73xl0" 2 

1.71xl0 -2 

1.62 xl0“ 2 

1.61 xlO" 2 

1.57xl0" 2 

12 C 

1.60xl(T 9 

1 58xl0 -9 

9.76x10"® 

9.62x10-® 

4.18x10-® 

16Q 

4.24x10"® 

4.19x10-® 

5.51x10"® 

5.45x10-® 

2.16xlO" 10 

20 Ne 

8.55xl0~ 12 

8.46xl0" 12 

7.05X10" 11 

6.96X10" 11 

2.92xlO" 10 

24 Mg 

4.38x10" 12 

4.33xl0" 12 

3.25xl0 -11 

3 21xl0 -11 

2.00x10"® 

28 Si 

1.92xl0 -11 

6 39xl0 -14 

1.92xlO- 10 

8.10xl0- 22 

2.95x10"® 

32 S ' 

4.94x I0" 11 

4 42xI0" J2 

5.34 xlO" 10 

2.14X10 -17 

2.08x10-® 

36 Ar 

2.60x 10 _1 ° 

1 26xlO -10 

2.46x10"® 

1.73xl0" 13 

1.59xl0" 7 

40 Ca 

1.19x10-® 

1 12x10-® 

1.90x10'® 

6.44x10-® 

4.18xl0 -7 

44 Ti 

4.73x10-® 

4.68x10“® 

1.71 xlO" 7 

1 69 xlO -7 

3.26x10-® 

48 Cr 

188x10"® 

1.51x10-® 

2.19x10-® 

1.26xl0~ 14 

2 55 xlO" 6 

52 Fe 

1.59xl0 -s 

1.61xl0" 5 

9.26xl0" 7 

1.35x10-® 

4.50x10-® 

56 Ni 

1.66x 10 -2 

1.66xl0- 2 

1.67xl0 -2 

1.67X10" 2 

1.64 xlO' 2 

60 Zn 

6.50x10-® 

6 39x10-® 

4.23x10-® 

4.15xl0" 5 

2.97xl0 -4 


larger nuclear networks, where the potential for improvement in speed and size is even 
greater. For larger networks, reduction m the number of nuclei by a factor of 2-4 could 
result in a increases in network speed of a factor of 5-10 because of the nonlinear relation 
between matrix size and the length of time to solve a matrix equation. 
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